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Gluon propagator and zero-momentum modes on the lattice * 
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We investigate the propagators of 4D SU(2) gauge theory in Landau gauge by Monte Carlo simulations. To be 
able to compare with perturbative calculations we use large /9 values. There the breaking of the Z(2) symmetry 
causes large effects for all four lattice directions and doing the analysis in the appropriate state gets important. 
We find that the gluon propagator in the weak-coupling limit is strongly affected by zero-momentum modes. 

1. INTRODUCTION 

Starting with Ref. there has been a num- 
ber of nonperturbative lattice studies of the gluon 
propagator in Landau gauge. However, the im- 
pact of zero-momentum modes on the propaga- 
tors has not been analyzed. Recently one of 
us has shown that zero-momentum modes may 
strongly affect gauge-dependent correlators 
Therefore, we have performed simulations in 4D 
SU2 lattice gauge theory to clarify this issue. 

To be able to compare quantitatively with per- 
turbative calculations large P values must be 
used. There the propagators, which gauge fixing 
effectively makes very nonlocal objects, become 
sensitive to the broken Z(2) symmetry states of 
the deconfinement region. For comparison with 
perturbative results the appropriate one of these 
states has to be selected. 

We use the Wilson action, periodic boundary 
conditions, fields Of^{x) — ^{Ufj^x — Uj^^) and the 
Landau gauge. The propagators considered are 
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rM(P, ^) = E Tr(OM(P, t + T)d^{^p, t)) (1.1) 



where 0^(p,r) = e'P'^C>(^> Choosing 

p — (0, OjPa) where p3 = j^p the transverse prop- 
agator is defined by Fx = hiTi +^2)- 
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Figure 1. Time history of Polyakov loops and 
(r — 0) propagators at /3 = 10 on 4"^ x 8 lattice. 

2. BROKEN Z(2) SYMMETRY 

In our simulations on lattices of sizes 4^ x 8, 
8"^ X 16 and 16^ x 32 we have used Polyakov loops 
= T^x^x^ Tr n^,^, U^^xto monitor the break- 
ing of the Z(2) symmetry. According to the two 
possibilities for each direction (i.e., > or 
Pf, < 0) there are 16 states : (+ + ++), (+ + H — ), 

The time histories of the Pf^, Tt and r4 in 
Fig. 1 illustrate that the are, in fact, strongly 
affected by the indicated states. This phe- 
nomenon has been observed on all lattices consid- 
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Figure 2. rT(j5, r) in various states for p — I, 
(3^10 and 163 x 32 lattice (solid line is (O)). 



ered. From the numerical analysis we find in more 
detail that the observables take different values in 
different states. As one can see in Fig. 2, (++++) 
gives results consistent with perturbation theory 
while otherwise large deviations occur. 

3. FORMS OF PROPAGATORS 

In the limit of large f3 = one obtains in 

lowest-order approximation 



2T/4^4sin^^+4sin^2^ 



(3.1) 



provided that 7^ 0. To handle we split off 
the zero-four- momentum part C^j, = ^ ^m(^) 
so that the fields decompose as (0, t) = + 
50f_t{T) and becomes 



r^(o,T) 



Tr(C2) + R^{t) ; 



(3.2) 



Mr) = ^Y.^i-{Sd^{t + T)Sd^{t)) (3.3) 
^4 t 



The evaluation of ( p.3[ ) by collective-coordinate 
methods only works in the lowest-order approxi- 
mation and leads to 



r,(0,r)^Tr(C^) + |f^^ 



2Vi ^ 4sin2f 



4. ZERO-MOMENTUM MODES 



(3.4) 



Now we restrict the considerations to data of 
the (+ + +-|-)-state. First we note important 
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Figure 3. Time histories of C^, 501, ^^(1) ^^'^ 
X'^{2) for /? = 10 and 16^ x 32 lattice. 



properties of by comparing the time histories 
of C^, 50 ^{t) and X^{p) = Re ^ e'^^-^O^ix) 
with p = (0, 0, -j^p) and p^O, i.e. of an example 
of a nonzero-momentum part of the fields. Fig. 3 
shows behaviors of components C° , SO"^ and 

(where = Ea('^V2)C^ etc.). For JO^, X-{1) 
and X^{2) uniform Monte Carlo noise is seen. 
Its magnitude is larger if the three-momentum 
involved in the particular quantity gets smaller. 
For C° in addition to such noise, with magnitude 
comparable to that of surprisingly large 

variations are observed. These variations exhibit 
different patterns for different p, while for differ- 
ent a we find essentially the same pattern. 

In the decomposition of the fields the parts 
and SOfj, obviously behave quite differently. 
The observed large variations of appear to be 
a characteristic consequence of zero-momentum 
modes. 

In addition to Tt{0,t) and Rt{t) also the 
quan tity Tr(C|n) needed for the comparison with 
(|3.4| ) has been determined in the simulations. 
From Fig. 4 it is obvious that our numerical re- 
sults for the transverse propagator with p — Q 
agree reasonably well with that of our lowest- 
order calculation (3^). It is seen that the zero- 
momentum part Tr(C|n) is large as compared 
to the rest, which is related to the large varia- 
tions of mentioned above. It is also apparent 
that while Tt{Q,t) exhibits unusually large er- 
rors, which stem from the large variations of C^, 
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Figure 4. rT(0, r) and Rt{t) compared with 
( ^^ ) (curves) and Tr(C^) (constant line in upper 
Figure) for /3 = 10 and 8^ x 16 lattice. 



for Rt{t) one gets errors of usual size. Analo- 
gous observations as for the 8"^ x 16 lattice with 
/? = 10 have been made on 4"^ x 8 with (3—10 
and on 16^ x 32 with /5 = 10 and /3 = 99. Thus 
our numerical results confirm the description we 
have given. 

The importance of zero-momentum modes is 
quantified by the fact that Tr(C^) is large as 
compared to Rf^iO). To compare different lat- 
tice sizes and different values of f3 the respective 
values are to be multiplied by Va/g^ (which ob- 
viously cancels the extra factor implicit in our 
definitions). It turns out that the values of 
^ = (V3/5^)Tr(C|.) are much larger than 7 — 
(V3/5^)(i?T(0) - i?r(i4/2)). From Table 1 it is 
seen that for increasing /3 this feature gets even 
more pronounced. The same holds for increasing 
lattice size. Not only C gets larger but there is 
also an increase of the ratios C/7. 



Table 1 
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Figure 5. Misleading results from determinations 
of effective masses from eq. (4.1) for p = 0, /3 = 
10. Curves are based on eq. (3.4). 



To show that the usual determination of ef- 
fective masses to(t) can be misleading if zero- 
momentum modes are present we have also cal- 
culated m(T) from 



cosh(m(T)(T + 1 



rT(0,r+l) 



(4.1) 



cosh(TO(T)(r- ^)) rT(0,T) 

From Fig. 5 it is seen that even at very weak 
coupling, where the situation clearly is described 
by (3^), a gluon mass is imitated. Its decrease 
with lattice size actually reflects the increase 
of the zero-momentum contribution. Because 
m(r)L^ shows little dependence on lattice size, 
the present results may also be considered from 
the point of view of finite temperatures where 
screening masses are determined 

To summarize, at /3-values considered here 
such masses are fake. The question of the role 
of the zero-momentum modes at smaller values 
of (3 (in the physical region) needs further study. 
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